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1. INTRODUCTION 

The rechargeable batteries become widely spread nowadays due to availability in markets in 
different sizes, as well as can be used for industrial applications or people needs. Another factor contributes 
to its spread is the environmental issue and the low-cost of this energy. However, theirs long life cycle and 
high density of its energy make it dominate in many applications, such as power grid, electrical vehicle, 
which have been integrated as an energy storage system [1]. In recent years, the rechargeable battery 
improved by using developed technologies such as limit gas emission and enhanced in mass and energy 
density. It is well known that the battery converts the chemical energy to electrical energy during discharge 
and the opposite during the charging process, due to this reason, modeling the battery becomes a necessity to 
control and safeguard the systems. 

Modeling the battery can give important pieces of information about its internal state. The internal 
state represented by the remaining energy inside the battery, which is signified by the state of charge (SOC). 
The main objective of modeling is to monitoring the internal parameters of the battery, among them the 
voltage, current, and SOC [2, 3]. Therefore, an accurate battery model is vital to forecast the battery 
characteristics and can avoid the overcharge or the over-discharge. In literature, there are a lot of studies on 
battery modeling. Some of them focus on electrochemical modeling [4, 5] and other about equivalent circuit 
models or combined models [6, 7]. Since the equivalent circuit models have been based on electrical 
components to suit chemical behavior it is applied widely for simulation. Thevenin and second-order battery 
model are popularly used for modeling the battery, but before selecting the battery model, it is necessary to 
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follow this steps. First, settle an experiment for standard tests to estimate the battery parameters. The second 
step is selecting an accurate model with minimal errors. 

Matlab/Simulink is a powerful tool for simulation; it provides a lot of blocks and tools such as 
functions blocks, and also strong libraries, as SimPowerSystems library. However, the model is built based 
on the electrical equations and developed in Simulink in order to give accurate simulation results. The paper 
is investigated on the equivalent circuit model build in Matlab/Simulink. The second-order Thevenin battery 
model was designed and the parameters of the battery model determination from experimental data. 
Moreover, a proposed model can be easily linked to other circuit blocks for comparison between real 
estimation and simulated data. Then, an extended Kalman filter algorithm (EKF) [8-10] was developed to 
enhance SOC estimation. Therefore, the coulomb counting method [11-13] has been used as a reference for 
SOC estimation. The paper is organized as follows: Section two introduces the proposed model and detailed 
by equations and blocks. Test bench and parameters extraction represented in section three, followed by 
outline the EKF approach in section four. The results of simulation are summarized in section five. Finally, 
the paper concludes. 


2. PROPOSED BATTERY MODEL 

The commonly used equivalent circuit for describing the battery model is Thevenin, which is 
representing the model by electrical components. It consists of an ideal battery voltage source (OCV), a 
series resistance (Rs), and a parallel capacitance (C) and an over-voltage resistance (R). The Thevenin model 
also named the first-order model, second-order model or third-order model depending on the number of the 
branch (RC) consist. The RCs branches are used to describe the behavior or transient response of terminal 
voltage. Figure 1 shows the selected equivalent circuit model. This second-order Thevenin model is built in 
Simulink and standard tests were done into the battery to determine its parameters. 


Figure 1. Equivalent circuit model 


Where V, represent the terminal voltage of the battery. J, is the current flow for charging or 
discharging (positive while discharging and negative while charging). OCV is representing the open circuit 
voltage, C, and R, represent the rapid transient voltage denoted as V,;, C} and R, form the second RC branch. 
R, represents the series resistance. 


2.1. SOC calculation 

State of charge is the famous parameter reflect the quantity of energy remain in the battery. It can be 
affected by some conditions as the temperature and the history of the cell. When the cell is fully charged the 
SOC take the value of 100% or 1, in opposite, when it is fully empty or discharged get the value 0. Numerous 
methods assume to determinate the SOC as the coulomb counting. It is the simplest approach based directly 
to the definition of SOC, the precision of this method depending on the initial value of SOC and require re- 
calibration each cycle of charge/discharge, mathematically expressed in [14]: 


SOC = SOC, — fat (1) 
bat 


Where the SOC, is the initial value, Cpa is the nominal capacity given by the manufacture. By 
converting the (1) to Simulink model for simulation purpose, Figure 2 shows the structure of SOC 


calculation. 
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Figure 2. Calculation of SOC 


2.2. OCV calculation 

Open circuit voltage is the equilibrium voltage. It is an important parameter to be identified which is 
the linked directly to the SOC. The relationship OCV-SOC can be represented with polynomial equation. 
However, a mathematical function can be establishing that relationship as shown in Figure 3. 


aE fe 


soc G ol OCV 
Fen 


Figure 3. Calculation of OCV 


2.3. RC values 

Figure 4 shows the structure of the RC branches, the value of R_1 , C_1, R_2 and C_2 are tabulated 
in look-up table due to their relationships with SOC. The model determines the best value for parameters by 
using interpolation-extrapolation method. 


Switch4 


Figure 4. RCs values 


State of charge estimation by using extended Kalman filter based on improved open... (Souaihia Maamar) 


4 g ISSN: 2252-8792 


2.4. Voltage of RCs branches calculation 

The branches RCs are corresponding to the transient response of battery terminal voltage. In order to 
build a subsystem responsible to represent the RC parallel network, in Figure | appear the voltage V, and V3, 
which are corresponding to transient voltages | and 2, their equations can be expressed as: 


LAG eee L (2) 
dt Cy, RG 


die AD a (3) 
dt C2 Rolo 


Thus, the blocks used to realize the simulation of these equations can be arranged as illustrated in 
Figure 5, where the first and the second RC network are V, and V> respectively. 


Divide7 - 
Divides Vi x 
R1 
Divide5 
1 
; 
Divide3 


Divide4 V2. 


Figure 5. RCs branches voltages 


2.5. Vs calculation 

It is necessary to mention that the series resistance Rç represent the difference between the drop 
value of voltage after pulse discharge, and it also shows different values when transferring between charge 
and discharge operations. The value of the voltage across it can be calculated as e (4), the structure block in 
Simulink is illustrated in Figure 6. 


kekeh (4) 


Product sa 


Figure 6. Vs voltage calculation 
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Finally, the terminal voltage (V,) is calculating by forming (5) as denoted below and the structure 
final of the proposed model formed by subsystem blocks is illustrated in Figure 7. 


V, = OCV(z) -V, — V, — V, (5) 


RC look-up 


RC voltage 


Figure 7. Proposed simulation model in Simulink 


3. BATTERY TESTS 

As it is illustrated in Figure 8, the schematic blocks show the experiment test bench structure. The 
bench consists of a LifeP04 battery, which has a nominal voltage of 3.2 and a capacity of 10Ah. The battery 
is discharged by electronic load and the data transferred to the host computer using Arduino mega with 
Matlab/Simulink in a room at 25°C. 


Host Computer 


communicatio board 


Sensors and 


Figure 8. The schematic of battery test bench 


3.1. Series resistance 

The series resistance Rs can be identified from experimental data as illustrated in Figure 9. The 
internal resistance can be calculated in the drop voltage between U1 and U2 during relaxation period (I=0A) 
by using (6). The obtained values tabulated in look-up table as shown in Figure 7. 


= U2-U1 (6) 


IL 


Rs 
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OCV at 90% 


OCV at 80% 


Voltage(V) 


2000 2500 3000 3500 4000 
Time(S) 


Figure 9. Transient response of terminal voltage 


3.2. Calculation of R1, R2, C1 and C2 

The values of R;, R2, C; and C, are determined from the transient voltage between U2 and U3. The 
least square method algorithm through curve fitting Matlab tool has been used to compute each parameter of 
RCs branch network. The obtained values are tabulated in look-up table as mentioned in Figure 4. The 
equation used to determine these values is denoted as: 


=¢ 


Uz = l, * R, * (e(raex)) + I, * Ry * (eraa) (7) 


3.3. OCV relationship 

The open circuit voltage is the terminal voltage of the battery after pulse discharge applied where the 
voltage reaches the equilibrium state. The values of OCV linked to SOC. The relationship between them is 
represented by polynomial function as shown in Figure 10. It is constructed by different experimental points 
fitted with polynomial function. The fitted polynomial function is defined as: 


OCV(z) = a,z° + apz* + azz? +4427 + asZ + A, (8) 


where z denote the state of charge (SOC) and the parameters a, to a, are tabulated in Table 1: 


Table 1. Values of the polynomial parameters 


Parameters Values 
ay -6.704 
az 16.72 
a3 -14.01 
ag 4.69 
as 0.03277 
a6 3.367 


Experimental data 
® polynomial fitted data 


Open circuit voltage (V) 


soc 


Figure 10. The OCV-SOC relationship 
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4. EXTENDED KALMAN FILTER APPROACH 

The Kalman filter is the most technology known on tracking systems, used also for the estimation, 
there exists for linear systems and has adopted for nonlinear systems under the name extended Kalman filter 
(EKF) [15-19], which is utilising for prediction into dynamic fields, for example: tracking, navigation, and 
also in global position. So recently has used in battery field. The Kalman filter consists of mathematical 
recursive equations, which allows the repetition of the calculi each step, the output of the system depend 
directly to the previous input, so the state vector relies on the present output to its input. Because of the 
battery behaviour and its dynamic. It can be considered as a nonlinear system for this part, due to the voltage 
and current sensors from another part. In this system, the output is the voltage, and the input is the current, 
also an observer is designed to calculate and evaluate the state each sampling time. The observer is designed 
to track each parameter; every parameter is representing into the system by a specific state variable. The 
system is based on the equations of the first order RC equivalent model, which cited in battery model chapter. 
The EKF technique is applied to suit the behaviour of the battery and extract the parameters each iterations. 
By supposing the state variables for the battery are: 


x=[V, V, SOC]? (9) 
with T is the transpose matrix. 


To provide a background, it is necessary to briefly describe the extended Kalman filter. For the state 
vector has been used a nonlinear continuous equation, and for the observation vector used a nonlinear 
discrete equation. The input is chosen as u(t)=IL and the output is y(t)=Vt, the equations of the system and its 
observer become [20]: 


Raed oe 
y=g(x,w) +w (10) 


The terms w and v are random distributions Gaussian representing the error caused by the sensors 
and measurements. In the previous functions f (x,u) and g (x,u) linearization is needed, so each sampling step 
is carried by the Jacobin function as follow [21]: 


T 


e R220) 0 
= 7 _ [44 aocv 
an = 0 e Pict of Ce = | aa, al 
0 0 1 
penis (11) 
R,(1—e F262) 
T 
Be =|R,(1—e 11) |> De = [Ro] 
ai 
Cn 
The EKF algorithm follows the steps: 
a. Prediction: 
Ky = Ar * yy + Bp * Tha (12) 
Pg = AkPk-14 +R (13) 
b. Measurement update 
K, = P, C™[CP, CT + Q]7? (14) 
Ry = Rg + KV. CK) — HCH) (15) 
P, = [I — KC] P7 (16) 


where Q and R are the variances of the noise v(k) and w(k) respectively, and 7 is the identity matrix, £p is the 
update state estimated from the previous one. 
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5. SIMULATION RESULTS 

In order to validate the proposed model, the simulation results of continuous discharge and pulse 
discharge tests are used to compare with the corresponding experiment results. Therefore, the structure of the 
schematic made for simulation is shown in Figure 11. The EKF is integrated to estimate the SOC. 


batt_time batt_current] 


Current 
Initial_soc 
Ah Method 
Terminal voltage 
Vt vs Vt_simulated 
E 
Add Error EK Function 
Parameters Memory 
1 i | 
batt_time batt_volt} J 
| : 
Voltage 


Figure 11. Schematic blocks of EKF SOC estimation 


By using the experimental data, such as current and voltage, where the input is the current-time 
vector and the voltage-time as an output of the model. As shown in Figure 12, the voltage curve of simulation 
and experiment results are well tracked. The RMS error of voltage is 25mv, which is 2% to the terminal 
voltage of the battery. This shows good agreement between the model and experiment, so the model can 
perform accurately. Similarly, the error of estimation the voltage is shown in Figure 13. 
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Figure 12. Comparison between simulated and experiment Figure 13. The error between simulated and experiment 
data of pulse discharge voltage 
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The improved battery model is suitable for SOC estimation, the hybrid pulse power test was 
implemented to figure out the accuracy of this model. The parameters of the battery are updated via look-up 
table and extrapolation, Figure 14 presents the comparison of real data with model simulation data. The 
figure shows a good match between real and simulation in SOC. As shown in Figure 15, the simulated SOC 
match the real SOC while the discharge current applied. The EKF SOC estimator instantly diminishes the 
error, while error decreases below 3% in continuous discharge. In our case, while the correction takes a few 
milliseconds. 


Estimated SOC Estimated SOC 
1 Real SOC 7 Real SOC 


O 0.6 
0.5F 
0.4; 
0.3 
0.2 
a 05 1 15 2 2.5 0 0.5 1 1.5 2 2.5 3 
Time (S) x10 Time (S) x10 
Figure 14. Estimation of SOC of pulse discharge test Figure 15. Estimation of SOC of pulse discharge test 


6. CONCLUSION 

This paper has been proposed as an accurate battery model to characterize the dynamic behavior of 
battery based on Matlab/Simulink. The EKF algorithm is applied to optimally estimate the SOC of the cell, 
the battery model has been based on the second-order Thevenin model, where has been implemented and 
explained in details. This model has been developed and the accuracy has been proven with experimental 
results. It is applicable for other battery models, while it is simple and not involves complex computation, as 
well as, the response of EKF algorithm in order of a few seconds. 
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